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Abstract 

In this paper I first address the question of whether the seat of the power radiated by an 
antenna made of conducting members is distributed over the "arms" of the antenna according 
to —J • E, where J is the specified current density and E is the electric field produced by that 
source. Poynting's theorem permits only a global identification of the total input power, usually 
from a localized generator, with the total power radiated to infinity, not a local correspondence 
of —J • E d 3 x with some specific radiated power, r 2 S • f dQ. I then describe a model antenna 
consisting of two perfectly conducting hemispheres of radius a separated by a small equatorial 
gap across which occurs the driving oscillatory electric field. The fields and surface current are 
determined by solution of the boundary value problem. In contrast to the first approach (not a 
boundary value problem), the tangential electric field vanishes on the metallic surface. There is no 
radial Poynting vector at the surface. Numerical examples are shown to illustrate how the energy 
flows from the input region of the gap and is guided near the antenna by its "arms" until it is 
launched at larger r/a into the radiation pattern determined by the value of ka. 
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I. INTRODUCTION 



This paper is a didactic discussion of how the electromagnetic energy radiated by an an- 
tenna emerges from a localized source, is guided by the antenna's conductors, and ultimately 
shakes free to form the radiation described by the asymptotic Poynting vector. How this 
happens and what is the true seat of the power radiated is understood by many, but not by 
all. I hope that an analysis of the rights and wrongs and the treatment of specific examples 
prove useful. 

The question of how electromagnetic energy is transported along a system of conductors 
dates back to J. H. Poynting's 1884 paper, On the transfer of energy in the electromagnetic 
field, in which he enunciates his theorem and discusses various examples. 1 Poynting does not 
consider radiating systems, but he was very clear on the transport of the electromagnetic 
energy associated with current-carrying wires and quasi-static circuits such as a discharging 
capacitor. 

In modern notation Poynting's theorem takes the form, 

— / J • E d 3 x = / tt d 3 x + <£ n-S da (1) 

Jv Jv at Js 

Here S = E x H is the Poynting vector, E and H are the fields, J is the current density, 
and, with some qualifying caveats that need not concern us here, u = |(E • D + B • H) is 
the electromagnetic energy density within a chosen volume V bounded by the surface S. 
The physical interpretation of the theorem is that the left-hand integral represents the rate 
at which the given sources (J) supply energy to the electromagnetic fields. The volume 
integral on the right is the rate of increase of electromagnetic energy within V, while the 
surface integral is the rate at which energy escapes from V through the surface S. The 
theorem is basically a statement of conservation of energy. 

In the late 19th and early 20th centuries, radiated systems were investigated chiefly by 
specification of simple oscillating current and charge distributions. The asymptotic fields 
were found and the Poynting vector evaluated to give the radiation pattern. Performing the 
surface integral on the right in gave the total radiated power P in terms of the parameters 
of the source, e.g., its dipole moment. For an antenna with a given input current /, the result 
was often expressed in terms of a radiation resistance R ra d through the relation, P = I 2 R ra d. 

In 1922, Brillouin 2 pointed out an alternative method for calculating R ra d- For steady- 
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state sinusoidal oscillation of the source, the volume integral on the right in (P) vanishes for 
a time average; the input power is equal to the power radiated. The radiation resistance can 
be computed equally well by evaluating the integral over the sources on the left-hand side 
of (Q). This approach was taken up by many^ because of its appeal: the current of each 
element along the arms of an antenna was apparently the real source of the radiated energy. 

At that time (and to some extent today) many believed that energy was transported 
within the conductors carrying current. Poynting's statements to the contrary had been 
forgotten. The focus on specification of distributed current sources is widespread in many 
applications. For antennas the approach works, more or less. Often antennas have thin 
straight conducting elements. If such elements have infinitesimal cross sections, it can be 
shown that an oscillating current is distributed sinusoidally along its length. The assumption 
of sinusoidal behavior proved to yield reasonably accurate results for angular distributions 
of radiation from antenna arrays. Specifying a plausible current source J became standard 
practice. In such discussions no attention is paid to the fact that a radiating antenna is an 
electromagnetic boundary value problem in which the current distributed on the arms of the 
antenna emerges as part of the solution, not as input. Rather, an ethereal current distri- 
bution is postulated and the resulting fields calculated. This approach may be appropriate 
for electrons orbiting in a synchrotron or undergoing a quantum-mechanical transition in an 
atom, but not for a careful treatment of a realistic antenna with conducting surfaces. 

The antenna as a boundary value problem has, of course, an honorable history beginning 
in 1897 with Pocklington 5 who discussed the lowest mode of a perfectly conducting ring of 
wire (damped by radiation loss). Lord Rayleigh extended the work to the higher modes in 
1912, with special attention to the damping.-- The history is summarized by R. W. P. King 7 
in his Introduction. Noteworthy are the works of L. V. King 8 and Hallen 9 in the 1930s on a 
long cylindrical antenna with a solution by means of an integral equation. In the early 1940s, 
Schelkunoff, in a paper^ and a bookii, discussed perfectly conducting biconical antennas 
in which the cones act as wave guides. A spherical antenna with a planar circular gap can 
be viewed as a limiting case. 11 At the same time, Stratton and Chu analyzed three models 
of antennas: the circular cylinder— ,the sphere^, and the prolate spheroid.- 14 Subsequently, 
R. W. P. King and colleagues made extensive use of Hallen's approach for linear antennas. 
Their results and literature are thoroughly documented in King's book. 7 A sampling of 
the accumulating literature on cylindrical, biconical, and spherical antennas in the applied 
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journals is contained in references . 15 1 16 i 17 i 18 

In none of these works is attention paid to the fields and energy flow in the immediate 
neighborhood of the antenna. F. M. Landstorfer and co-workers are one group that has 
explored numerically field lines and power flow in the neighborhood of guiding conducting 
surfaces: a half-dipole antenna on a conducting plane,— wave-guide discontinuities and 
dipole antennas,— two-dimensional diffraction,— and examples of gain-optimized antennas.— 
Surprising eddies of power flow can be seen in some examples. 

In Section II, I discuss a thin center-fed linear antenna by specification of a sinusoidal 
axial current density of vanishing cross section, with no attempt to solve a boundary 
value problem. The two different approaches to determining the total power radiated are 
described. The question of whether the element of "input" power —J ■ E at a particular 
point along the antenna can associated with a particular segment of the radiated power 
is answered. I also make some comments here on reasons for the relative success of this 
naive approach. In Section III, I describe a particular antenna boundary value problem, a 
perfectly conducting sphere with a small gap at the equator, across which is an azimuthally 
uniform electric field. As already stated, this problem has been addressed by Stratton 
and Chu^ and by Schelkunoff 23 . These authors were interested in the impedance of the 
antenna, including radiation, but did not discuss the Poynting vector in the neighborhood 
of the antenna or its evolution to the far fields. Section III sets up the formalism; Section 
IV gives examples of the energy flow, the current on the surface, and the modification of the 
energy flow for non-vanishing surface resistivity. Section V contains concluding remarks. 
An appendix gives the expansion in associated Legendre functions of the electric field at 
the gap, a necessary ingredient for the multipole expansion of Section II. 

II. EXAMPLE OF LINEAR ANTENNA 

As discussed in the Introduction, a common approach to the emission of radiation by 
an antenna is to postulate (with greater or lesser justification) the sources as given current 
and charge distributions. The potentials and the fields are then evaluated in the usual way 
as integrals over the source with the appropriate Green function. We discuss briefly the 
very thin center-fed linear antenna in order to illustrate the confusion that can occur in 
evaluating the radiated power in different ways via Poynting's theorem and in arguing the 
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source of the power. 

The time dependence is assumed to be exp(-iut) with k = uj/c. Complex notation is 
used, with physical quantities as the real parts of complex expressions. The antenna is 
of length 2a, located on the z-axis on the interval —a < z < a. The current and charge 
densities on that interval are assumed sinusoidal and are 

J(x, y,z) — z Io sin(/ca — k\z\) S(x)S(y) 

p(x, y, z) = i— e(z) cos(ka — k\z\) S(x)5(y) 
c 

With the time dependence suppressed, the vector and scalar potentials in the Lorenz gauge 
at a point (x, y, z) are 



A = z- 



An 



T dz' sm(ka-k\z'\) e lkR /R 

J —a 

- r dz' e(z') cos(ka - k\z'\) e ikR /R 

J ~a 



Here Z = ^Jp /e is the impedance of free space, R = yj p 2 + (z — z') 2 , p = \Jx 2 + y 2 , and 
e(x) = ±1 for x > 0, x < 0. To evaluate the scalar potential's contribution to the electric 
field we need 

-V(e lkR /R) = V'(e ikR /R) 

The electric field is the sum of the negative time derivative of the vector potential (Ei) and 
the negative gradient of the scalar potential (E 2 ) : 



E x 



^ZqJo f a dz ,^ k gin(A;a _ k ^ykRj R 

4:71 J -a 



Eo = % 



Att 



f a dz' e(z') cos{ka - k\z'\) V'{e lkR /R) 



Taking only the ^-component of the gradient, we integrate by parts in E 2 to obtain 



R 



2z 



■ ZqIq 
An 



f . a p ikR a 

{e ikri / ri + e lkr2 /r 2 ) - / dz' -^-[e(z') cos(ka - k\z' 

J -a R dz 1 



(2) 
(3) 

(4) 



Here T\ = yj p 2 + (z — a) 2 and r 2 = \J p 2 + [z + a) 2 . The derivative of the charge density has 
a delta function contribution at the origin from the derivative of e(z') so that, when we add 
in the transverse contribution we have 



Eo 



-Ex 



+ i 



■ Z I 



An 



(e ifcr 7n + e lfcr2 /r 2 - 2 cos(A;a)e jfc 7r) 



+ r dz' e(z') cos(ka-k\z'\)V ± (e lkR /R) 

J— a 



(5) 



Here r = \/ p 1 + z 1 is the distance from the origin. The total electric field is therefore 

. Z I 



E 



47T 



z (e ikri /n + e ikr2 /r 2 - 2 cos(A;a)e ifc 7r) 



+ / dz' e{z') cos{ka - k\z'\)V ± {e ikR /R) 



(6) 



If we wish to evaluate the time-averaged power input Pi to the electromagnetic field, 



Pa = - J d 3 x (J • E) 



we need only the z component of the electric field on the axis. Explicitly, we have 

P 1 = -~Re f dz I z E* z (p = 0,z) 

2 J —a 

With 



(7) 



(8) 



E z (0,z) 



47T 



e ik\"-z\/\ a _ z \ 



+ e 



ik\a+z\ 



\a + z\-2 cos(ka)e iklzl /\z 



(9) 



we find 



Pi 



z i 2 r 



in j- 

2a 



a 2 — z 



f 

/ dz sin(fca — k\z\) 

J —a 

L sin(fca) cos(kz) — (ajz) cos(fca) sin(kz) ] 



(10) 



The first line in the integrand is recognizable as the current; the second is the electric field 
after the use of some trigonometric identities. 



By Poynting's theorem the input power Pi is equal to the integral of the outward normal 
component of the time-averaged Poynting's vector S = Re(E x H*)/2 through any closed 
surface surrounding the antenna, in particular a sphere of large radius centered on the 
antenna. The well-known resul t 24 ' 2 ^ for the Poynting- vector power P 2 is 



Z Iq r 1 | cos(ka cos 9) — cos(/ca)| 2 

F2 = ^rJ-r { cos0 > i - cos 2 e 

We write with malice aforethought ( = a cos 9 and transform P 2 Into 

cos(/c£) — cos(ka)\ 2 



Po 



Z P a r- 



(11) 



(12) 



4tt J- a "" a 2 -C 2 
Comparison of Pi and P 2 suggests several questions: 
1. Are the two expressions actually equal? The answer is yes. Each can be transformed into 



the same uninformative sum of constants, logarithms, and sine and cosine integrals 24 that 
undulates SIS db function of ka. 

2. Does the variable z in Pi correspond in any way to ( = a cos 8 in P 2 and so connect 
the increment of power at position z on the antenna to the increment of radiated power at 
angle 9 ? The answer is no for several reasons. First of all, the integrand in Pi is a spatial 
scalar, without direction, whereas the integrand in P2, although a scalar, is the dot product 
of the Poynting vector with a unit radial vector whose direction changes with angle. 

Secondly, a comparison of integrands shows smaller and larger differences between the 
two "angular" distributions depending on the value of ka. As a matter of fact, for ka = 
(2n + l)vr/2, n = 0, 1, 2, . . . , the integrand of (jlUj) as a function of z/a is equal to that 
of (jllj) for cos 6 = z/a, but for other values of ka the integrands can be vastly different, as 
shown for example in Figure ^ for ka = 2n. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

z, cos 



FIG. 1: Normalized integrands for the two forms Pi and P2 for the total radiated power from a 
linear antenna with ka = 2tt. The solid curve is the angular distribution of radiation in cos 9 from 
the Poynting vector expression P 2 . The dotted curve is the integrand J • E of Pi versus z/a. 

3. Can the expression for P\ justify the claim that the power radiated by an antenna has 
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its origins at various points along the length of the antenna? The short answer is no. 

First of all, the fields at any point in space receive coherent contributions from all elements 
of the current along the whole antenna. Kingi (p. 562-5) discusses a center-fed antenna with 
ka = 3tt/2, imagining it as the limit of a line of closely spaced individual generators that 
produce the sinusoidal current distribution. He shows that in a certain sense one can identify 
each of the three separate "ears" of the angular distribution at infinity with a corresponding 
half-wavelength segment of the current distribution. But he is at great pains to make clear 
it does not mean that the fields in each separate angular region are "maintained entirely 
and exclusively by the current in [the corresponding] third of the antenna. . . . The field at 
any point in space is maintained by all currents in the entire antenna." (op. cit., p. 565). 

Secondly, as stressed in the Introduction, a source distribution (J, p) in the Maxwell 
equations does not define appropriately an antenna boundary-value problem. Antennas are 
almost universally constructed of conductors, usually very good conductors. An excellent 
approximation for the boundary-value problem is to assume perfect conductivity. Then 
there can be a surface current density, but no tangential electric field at the surface. There 
is therefore no J • E to integrate over the surface to find the power. You might say, well, in 
the real world conductivities are not infinite. Suppose we assume excellent, but not perfect, 
conductivity. In that situation things go the opposite way: There is now a tangential 
electric field at the surface, but the Poynting vector points inward! Energy flows inward 
into, not out from, the conductor; the resistance of the metal converts the electromagnetic 
energy into heat; the antenna robs power from the sourc e 11 ^ . 

As already said in the Introduction, most antenna experts seem to accept the idea that 
the antenna structure does not itself radiate (in the sense that segments of the conductor 
are the local sources of power), but rather guides the energy from the input source and 
launches it into its final radiation pattern. Why then is the approach of specification from 
the beginning of the current on an antenna so widely accepted and employed? One reason 
is that the most common antenna is linear, or an array of linear elements. For a thin linear 
antenna the zeroth order current is sinusoidal in amplitude and constant in phase with 
position. The solution 9 of the boundary value problem yields corrections to the current as 
a series in powers of where Q = 2ln(2h/r), with 2h the length of the antenna and r its 
radius. Values of Q ~ 10 ± 3 is a plausible range: a 100 MHz A/2 antenna with 2h — 1.5 m 
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and a diameter of 0.5 (1.0) cm has Q = 13(10). 

Calculations including first order corrections show that for short antennas (kh < 2) the 
real part of the current either dominates the imaginary part over almost the whole length 
of the antenna, or the relative phase changes little over the length. In either circumstance 
the radiation pattern of an antenna or antenna array will differ little from that computed 
with the naive zeroth order current. Even for larger kh, for which the zeroth order current 
changes sign along the length, the imaginary part of the current is generally small compared 
to the real part (except in the narrow regions where the real part vanishes). Examples of 
the current components on center-fed linear antennas of different kh are given by Hallen^ 
and King.— 

If reliable values of input impedance are desired, second order corrections to the current 
at the input point may be required.— But for radiation patterns of linear antennas the naive 
assumption of a real sinusoidal current is adequate in most circumstances. 

We now turn to a solvable model of an antenna to show in detail how the energy flows 
in its neighborhood and beyond. 

III. CENTER-FED SPHERICAL ANTENNA 

To illustrate the near fields and power flow in the immediate neighborhood of an antenna, 
we consider a spherical antenna of radius a with excitation provided by an electric field 
across a narrow equatorial gap between the perfectly conducting near-hemispheres. As 
noted earlier, this problem has been treated in some respects by Stratton and Chu^ and 
Schelkunoff. 2 - We follow the notation of Section 9.7 ff in my text.— The power source does 
not concern us except that it manifests itself as an electric field across the gap (and the 
accompanying surface current on the surface of the sphere). The source can be imagined in 
Schelkunoff 's language of conical antennas. Our antenna can be thought of as a limiting 
case of a biconical antenna with spherical segments as caps and an opening angle of almost 
90 degrees. The sides of the cones meet at the origin inside the sphere. It is there that an 
infinitesimal generator resides. It causes an electric field between the upper and lower cones 
to emerge in an azimuthally symmetric way as Eg across the gap. 

A. Fields and Surface Current 
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Outside the sphere the electric and magnetic fields can be described by multipole fields.— 
The magnetic multipole fields have a non-vanishing radial component of the magnetic field, 
while the electric multipole fields have a non-vanishing radial component of the electric 
field. Because the normal component of magnetic field vanishes at the surface of a perfect 
conductor, while the normal component of the electric field does not, we conclude that the 
spherical antenna has only electric multipole fields. Furthermore, because of the assumed 
azimuthal symmetry, only (£, m) multipoles with m = occur. The electric and magnetic 
fields for r > a that represent outgoing waves at infinity are 



E 



H 



-zoi: 

£ 



a{£) 



a(£) 



D e (kr)-^Y e , 



(13) 
(14) 



Here k = uj/c, a(£) is the £ th electric multipole amplitude, h!p is the spherical Hankel 
function of the first kind, and the function DAx) is 



2M*) = ~[s 



h 



X) 



X 



+ 



dh { p{x) 
dx 



(15) 



With dY i>0 /de 
written as 



(21 + l)/47r P}( cos 8), the 9 component of the electric field can be 

a(£) D e (kr) P}( cos 9) 



E, 



2£+l 



(16) 



, y An£(£ + 1) 

To determine the multipole coefficients we must equate this component of the field at 
r = a to the expansion found in Appendix A for the field across a symmetric gap defined 
by angles ±e on either side of 6 = ir/2. Equating (jTUjl with (jA7|) yields 



a(t)D t (ka) 



V 



7i(2£ + 1) P t ( sine) 



(£ odd) 



(17) 



Z a\ £(£+!) e 

For a symmetric equatorial gap, the even £ multipole moments vanish. The surface current on 
the hemi-spherical conductors is given by K = f x H(r = a+). There is only an azimuthal 
component of magnetic field (|14p. The surface current is therefore in the 9 direction: 



K e = -HJr = a+) 



£ odd 



\\Sl(fA) aW ^(H^oostf) 



(18) 
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B. Total Power Input from the Source 



The oscillating electric field and the associated magnetic field in the gap produce a radial 
power flow at r = a. This time-averaged input power is given by the integral over the 
segment of the sphere occupied by the gap of Re(E x H*)/2: 



1 / r 27T 
^Pinput = -Re [ a 2 d4 

l \ JO J — sin e 



d{ cos 9) E Q m 



(19) 



In Appendix A it is shown that the field Eg in the gap region (— sine < cos# < sine) is 

V 

E e (r = a) 



2ae sin 9 

With HI taken from (fTHj) . the input power is therefore 



Pinput 



/sin e 
d{ cos 9) Re 
- sin e 



i V 



2ae sin^ 



E 



odd 



(2£ + l) 



a*(£) Hf\ka) P}{ cos 9) 



\ Att£(£ + 1) 

With the definition of the associated Legendre function, 



(20) 



p/( cos 9) = - sin^ dP e ( cos9)/d( cos 9) 



the integral is elementary. The result for the input power is 



^input = Re I ' 



irVa 



E 

l odd 



^4§7TT) a * W ^ (sine) h?){ka) 



(21) 



Comparison of the terms in this sum with a(£) given by (|17|) shows that 

Z a 2 



Pi 



input 



W)\ 2 Re[-iD t (ka) hf\ka) ] 



(22) 



l odd 



Now the Wronskians of the spherical Bessel functions can be used to show that 

-iDAx) h? ) (x) 



i' 2 dx x 

Therefore the time-averaged input power becomes the standard multipole expression, 

Z 



(23) 



33 



P 



input 



2 k 2 



E 



(24) 



l odd 
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Note, however, that this result for the power is found as the input power at the source, not 
from integration of r 2 times the asymptotic radial Poynting's vector over all angles. 



IV. POYNTING VECTOR AND LOCAL ENERGY FLOW 



A. Poynting vector 



The fields around the spherical antenna are such that the Poynting vector has components 
in the radial and transverse (9) directions. In terms of the fields in Section III. A, the time- 
averaged Poynting vector is 

D Z 



Re r R{t,e,r) P}( cos0)P/( cos( 



" odd i,£' 

+ T(£,£',r) P £ ,( cos 0)P e ( cosi 



(25) 



The coefficients R(£, £', r) and T(£, £', r) are 



R(£,£',r) = — 
T(£, £', r) 



(2£+l)(2f + l) 



4tt^ ££'{£ + !){£' + I 
+ 1) h [ p(kr) 



a*{£')a{£) hf{kr)D t {kr 



R(£, £', r) 



(26) 
(27) 



kr Di(kr) 

Before discussing the flow of power near the antenna, we consider two limits of the real 
part of the coefficient R(£, £', r), when £' = £ and kr ^> £, £'. Using (j!5|) and the asymptotic 
forms of the spherical Hankel functions, we find 

1 (2£ + 1) 



kr > 



Re[ R{£,£,r)} 
lim Re[R(£,£',r)} 



47rPr 2 £(£ + 1) 

1 I (2£ + l)(2f +T)" 
Ank 2 r 2 \ ££'{£ + !)(£' + 1) 
i£(£+l) 



Re[(iy~ E a*(£')a 



kr 



R(£,£',r) 



(28) 
(29) 
(30) 



Note that (|2*Hj) holds for all r > a. Thus if the radial part of (|25J) is integrated over a sphere 
of radius r, the orthogonality of the associated Legendre functions, 

+ D 



P},{ cos 9)Pj{ cos 9) d( cosi 



2£+ 1 



-5 pi 
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plus (J28j) leads directly to (J23J, valid at any r > a. Obviously this must hold because of 
conservation of energy flow. 

The other limiting form shows that the asymptotic angular distribution of radiation 
depends on the real part of (iy'~ e a*(£')a(£) , whereas the radial power flow at nearer dis- 
tances can be expected to be different because of the more complicated structure of the full 
expression (|2bp. And at close distances there is power flow in the transverse (6) direction as 
well. However, the asymptotic form of T (j3*Uj) shows that that component of the Poynting 
vector falls off faster than r~ 2 . We explore these aspects immediately. 

B. Examples of energy flow near the antenna 

Figure 2 shows the Poynting vector pattern close to the antenna for ka = 1.0. Here, and 
everywhere else unless stated otherwise, the total angular gap is 2e = 7r/50. One sees that 
the power squirts out of the gap in a radial direction, with some transverse flow parallel to 
the surface, guided by the antenna surface. Figure 3 is an alternative display, with r/a = 10 
as well as r/a = 1.1. The appreciable transverse flow and sharply peaked radial flow near 
the antenna goes over into almost purely radial flow of the asymptotic angular distribution. 



\V7 
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FIG. 2: Distribution of Poynting vector for ka = 1.0 and r/a = 1.1. Bases of arrows are observation 
points, lengths of arrows are relative magnitudes. Close to antenna, power squirts out from gap 
with some transverse flow guided along surface 
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FIG. 3: Normalized distributions in cos# of radial and transverse power flow dP r /d{ cos#), 
dPt/d( cos 6) for ka = 1.0 and r/a = 1.1, 10. Close to antenna, the power is peaked strongly in the 
equatorial plane, with some transverse flow. At r/a = 10, the angular distribution is essentially 
the featureless asymptotic form, with negligible transverse flow. 

Figure 4 and r/a = 1.1, 2, 10 in Figure 5. Note that the behavior of the power flow close 
to the antenna is essentially independent of the value of ka, a consequence of the source 
being a small equatorial gap. The flow rearranges itself at moderate and large distances, 
however, into very different asymptotic forms for the different values of ka. The antenna 
surface guides the flow only fairly nearby. 

C. Examples of Surface Current 

The surface current density on the sphere is given by (|18|). Integration over the surface in 
azimuth for fixed cos 9 yields the distribution hheta = ^na sin 9 Kg, the total current flow 

14 



FIG. 4: Distribution of Poynting vector (times r 2 ) for ka = 3ir/2 and r/a = 1.1,4.0. Bases of 
arrows are observation points, lengths of arrows are relative magnitudes. Close to antenna, power 
squirts out from gap with some transverse flow along surface. At greater distance, flow is largely 
radial and begins to approach the asymptotic distribution. 

across the "latitude" circle at fixed 9, as the analog of the current I(z) for a linear antenna. 
[Actually, I sp here{z) = sin9 Ig is a more accurate analog.] The real and imaginary parts of 
Iq in units of V/Z are shown in Figure 6 for the example of ka = 1.0. The real part of the 
^-component of the current is not far from proportional to the current on a linear antenna 
with ka = 1, namely I sin(l — | cos#|). Note, however, the comparable imaginary part, in 
contrast to the thin linear antenna. 

The corresponding current for ka = 37r/2 is shown in Figure 7. Here the real part of 
the current displays some resemblance to the three half cycles of a sinusoidal current on a 
center-fed linear antenna of the same ka despite the differences in shape. 
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FIG. 5: Normalized distributions in cos 9 of radial (solid lines) and transverse (dashed lines) 
power flow dP r /d( cos#), dPt/d{ cos#) for ka = 3tt/2 and r/a = 1.1, 2, 10. Close to antenna, 
the radial power is peaked strongly in the equatorial plane, with considerable transverse flow. At 
r/a = 2, the equatorial peak has vanished, the transverse flow has diminished, and the radial flow 
begins to resemble the r/a = 10 angular distribution, which is essentially the asymptotic form, 
with negligible transverse flow. 

D. Modifications because of finite conductivity 

The treatment so far has been based on vanishing resistivity on the surface of the antenna. 
In such circumstances the only tangential electric field at r = a is in the equatorial gap. 
With small but non-vanishing resistivity, a small tangential electric field exists over the 
whole antenna. A perturbation approach can be used to find the fields and the inflowing 
Poynting vector at the surface.— In addition to the zeroth order electric field f|13j) . there is 
a first order tangential electric field at the surface, given by 




(31) 
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FIG. 6: Real (solid line) and imaginary (dashed line) parts of the total surface current Ig = 
2ira s'm8 Kg( cos#) in units of V/Zq versus cos# = z/a for an antenna with ka = 1.0. 

where R s = p/5 is the surface resistance, p being the resistivity and S being the skin depth. 
The magnetic field is given by (|14|) . In passing we note that finite resistivity causes a 
modification in the multipole amplitudes ()17|) : 

h { p(ka) 



a(£) -> a(£) 



Z D e (ka) 



(32) 



The time-averaged energy dissipation per unit area of the antenna can be calculated 
either from the real part of the inward Poynting vector at the surface, 

-r ■ - Re\E^ x H< >1 
2 L J 

or from one half the surface resistance times the square of the surface current density, 

dPloss 1 



-Rs\K e \ 



(33) 



dA 2 

Here \K&\ = , |, as shown in (fTK|l . Note that, in the Poynting vector expression, the 
zeroth order electric field does not appear; its contribution gives rise to the source-generated 
outward power flow at the gap. For orientation, we note that for copper at room temperature, 
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FIG. 7: Real (solid line) and imaginary (dashed line) parts of the total surface current Ig = 
2na sin 8 K$( cos 8) in units of V/Zq versus cos 8 = z/a for an antenna with ka = 3ir/2. 

R s /Z = 6.8 x 1CT 7 , 2.2 x 1CT 5 , at 1MHz, 1GHz, respectively. For aluminum alloys the 
numbers are 1.5 - 2.0 times larger. 

An example of the distribution of energy loss over the surface of the antenna is given in 
Figure 8 for ka = 3n/2. The dimensionless quantity displayed is the absolute square of the 
ratio of the surface current density at cos 9 to its value at cos 9 = sin e, the edge of the 
gap. This quantity is related to the power loss according to 

1 dP loss R s \K e { cosfl)| 2 

P rad d(cos9) 2irR rad | K e ( sin e) | 2 1 j 

The ratio of the total dissipative loss to the radiated power is found by integrating (J34|) 
over cos 9: 



Ploss 


Rs Y^odd 


„\a{e)\ 2 \ka h ( l\ka)\ 2 


p 

1 rad 




J2odd I a W 2 



(35) 



The factors \xh e (x)\ 2 in the numerator are polynomials in inverse powers of x 2 up to x~ 21 
(For I = 1 the polynomial is (1 + x~ 2 )). For ka << 1, the dominant multipole is £ = 1; 
the ratio of the sums is closely l/(A;a) 2 . An alternative way of looking at it is that, in the 
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FIG. 8: Distribution of resistive power loss per unit cos# for ka = 3ir/2. The plotted quantity is 
the absolute square of the ratio of surface current densities given in equation (|34[) . The losses are 
largest near the gap, as that is the region of largest current - see Figure 7. 

quasistatic limit, the dissipative losses are frequency independent (except for that in the 
surface resistance), while the power radiated is proportional to (ka) 2 . 

For the opposite limit, ka >> £, the spherical Hankel function factor can again be large 
compared to unity, causing successive terms in the numerator to fall off more slowly than 
those in the denominator. But since the initial terms in numerator and denominator are com- 
parable, the cumulative effect is not major and the ratio of sums is of order unity. Here are 
some examples: (ka; ratio) = (0.1; 119.6), (0.5; 5.62), (1.0; 2.21), (10; 2.62), (20; 2.66). 

V. CONCLUDING REMARKS 

The conventional way of specifying the plausible sources p, J of charge and current on 
antennas can yield reasonable radiation patterns, but fails to address the actual boundary 
value problems. By discussing first a conventional example of a center-fed linear antenna, 
I have, I hope, shown that the idea of a one-to-one correspondence between an increment 
of input power —J ■ Ed 3 x in or on the antenna and an increment of radiated power per 
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unit solid angle r 2 S • f is without merit. Only the total input power can be equated to 
the total radiated power. Treatment of a tractable, albeit stylized, perfectly conducting 
spherical antenna with gap excitation as a boundary value problem illustrates how the 
antenna's current distribution emerges as part of the solution. The fields right down to 
r = a permit calculation of the Poynting vector everywhere outside the antenna. Numerical 
examples show how the power flow originates at the gap, is guided near the antenna by the 
conducting surface of the antenna, and is launched toward its ultimate radiation pattern 
after a few multiples of r/a. With perfect conductivity there is no energy flow into or out 
from the antenna's surfaces. When resistive losses are included, a small radial component of 
the Poynting vector occurs at the surface of the antenna, directed inward into the conductor 
where it is dissipated in ohmic heating. 

The main message is that antennas are boundary value problems, that the "arms" of an 
antenna guide and launch the energy flow but are not its origin, and that initial specification 
of the current for a metallic antenna is not consistent (even if it may yield reasonable results 
in some situations). 

As a final note, the reader should know that there exists canned finite element analysis 
software called NEC4 for computation of fields and energy flow around antennas of relatively 
arbitrary shape. One option is a plot of energy flow lines in which the local density and 
direction of lines indicate the magnitude and direction of the flux. Some may prefer that 
presentation over my Figures 2 and 4. Personally, I rather like the action implied by the 
arrows in my figures. 
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APPENDIX A: SURFACE ELECTRIC FIELD AS EXPANSION IN ASSOCI- 
ATED LEGENDRE FUNCTIONS 



A spherical antenna has a gap in its perfectly conducting surface defined by cos#i < 
cos# < cos 02- The internal source of power creates an electric field Eg at r = a within the 
gap, uniformly in azimuth. Otherwise, Eg = on the surface. The multipole expansion (16) 
of Eg is in terms of the associated Legendre functions Pj. We thus require an expansion in 
those associated Legendre functions of the rectangular function, f(z) = [Q(z — Zi) — Q(z — 
z 2 )] where Q(x) is the Heaviside step function. We begin with the completeness relation 
on the interval (—1,1) in z = cos 9 for the Legendre polynomials Pe(z) and P}{z) = 



-VT^^dP t (z)/dz: 

1 00 

*(*-*0 = «E(«+1) Pt{z')Pi{z) (Al) 
We integrate (jAljl in z over the interval (— 1, z) to obtain 

1 00 rz 

Q(z-z') = + / P*(z)<fc (A3) 

For £ = the integral is 

\ £ PoW d2! = (1 + 2)/2 

and for t > 0, 



The final result can be found in Magnus, Oberhettinger, and Soni.— Substituting these two 
results into (|A3|) yields 

Q(z-z') = i±£ + > /r^gi^ti P ^)^( z ) (A4) 
We define the tangential field on the surface as 

E e (r = a, z, Zl , z 2 ) = [ G(z - z x ) - &{z - z 2 )\ (A5) 

V 1 — z 

where A will be chosen for convenience below. With the expansion (|A4|) for Q(z — z') the 
tangential electric field in the gap is given by 

00 (21 + 1) 

E e (r = a,z, Zl ,z 2 ) = A E ^ + 1} [ fl(*0 - Pfa) ] P}(z) (A6) 
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In our calculations we chose the gap to be relatively small and centered around cos 9 = 0. 
With z\ = cos(7r/2 + e) and zi = cos(7r/2 — e) we find, using (jA5|) . that the voltage V, 
defined as the integral of the electric field across the gap, is V = 2aeA. From the symmetry 
of the Legendre functions around 9 = tt/2, a symmetric equatorial gap implies only odd i 
terms in (jA6|) . The result for the tangential (^-component) electric field on the surface of 
the sphere and in the gap is 

V (2P + 1) 

E e (r = a, cosfl) = — £ jTj^Pii sine) P/( cos£) (A7) 

l odd \ ' 

In the limit of e — > 0, equation (jA7|) can be shown to be equal to V/a times flA2|) with 
z' = 0. 
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